7-hadron correlations as a tool to trace the flow of energy lost from hard partons in 

heavy-ion collisions 
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High transverse momentum (Pt) 7-hadron correlations are currently being regarded as the 'golden 
channel' for the study of the medium produced in ultrarelativistic heavy ion collisions by means of 
hard probes. This is due to several reasons, all linked to the fact that because of the smallness of 
the electromagnetic coupling a, the photon does not substantially interact with the medium and is 
expected to escape unmodified. Thus, a high Pt photon indicates a hard process in the collision 
independent of the position of the hard vertex. In contrast, there may not be a clear signal for a 
hard process involving strongly interacting partons if the production vertex is deep in the medium 
as both partons undergo substantial final state interaction. Equally important, if photon production 
by fragmentation can be separated experimentally, the photon provides almost full knowledge of the 
initial kinematics. In the present paper, these properties are used to demonstrate a distinguishing 
feature between two assumptions made in modelling the medium-modifications of strongly interact- 
ing high Pt processes: Loss of energy into the medium vs. medium modification of the partonic 
shower. Is it shown that 7-h correlations provide a very clean signature to distinguish the two 
scenarios. 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 

The main goal of studying medium-induced modifications 
to hard processes in heavy-ion (A- A) collisions is to do jet 
tomography on the medium [l|, 0, 0, El IE IE] • The under- 
lying idea is as follows: Since hard processes take place 
slightly before the creation of a soft bulk medium, high pt 
partons have to traverse the medium before hadroniza- 
tion, and interactions with the soft medium subsequently 
lead to an energy redistribution, resulting in a suppres- 
sion of observed hard hadron yield. The strength of the 
interaction with the medium reflects the density of the 
medium, thus one should be able to gain information 
about the medium density from the observed strength of 
the suppression. The chief observable considered so far 
has been the nuclear suppression factor of single hadrons 
Raa (cf. e.g. @, 0]) which is the measured yield in A- A 
collisions divided by the yield in proton-proton (p-p) col- 
lisions multiplied with the number of binary scatterings 
(i.e. the default expectation if there would not be soft 
processes forming a medium in an A-A collision). 

There are at the moment two main scenarios which 
are considered in models to explain the measured sup- 
pression: Energy loss and the medium-modified parton 
shower. In the energy loss picture, the approximation 
is made that a single parton carries most of the energy 
and momentum of the event. Note that this approxi- 
mation is justified for single inclusive hadron spectra or 
back-to-back correlations, but clearly not for fully re- 
constructed jets. In traversing the medium, this lead- 
ing parton loses energy to the medium and subsequently 
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fragments into a hadron shower in vacuum after emerg- 
ing from the medium with reduced ener gy. This pic- 
ture describes the available data well [gL lloL HH using 
e.g. the Armesto-Salgado- Wiedemann (ASW) implemen- 
tation [12j of energy loss. The energy loss model itself 
does not explicitly specify what happens to the lost en- 
ergy (it is computed in the form of gluon radiation, but 
clearly also radiated gluons would re-interact with the 
medium), but it is often assumed that lost energy excites 
Shockwaves in a hydrodynamical bulk medium, see e.g. 

[ilia eh, iia. 

In contrast, the medium- modified parton shower pic- 
ture is not restricted to the approximation that a sin- 
gle parton carries most of the energy and momentum, 
as it follows medium-modified evolution equations for 
the full sequence of branchings of the initially produced 
high virtuality parton. Only the nonperturbative part 
of fragmentation, i.e. hadronization, is assumed to take 
place outside the medium. There is a class of mod- 
els in which this evolution of a whole parton distribu- 
tion is followed analytically [IE 0, [2(| HI| and also 
Monte-Carlo (MC) descriptions of the in-medium shower 
[22I [23l . H3, [5E HE] based on the vacuum shower evolu- 
tion routines in PYTHIA [27| or HERWIG [H]. While 
these models can account for Raa well, usually the de- 
tailed comparison with more differential observables is 
not quite as advanced yet due to the greater complexity 
of these models. The important difference to the en- 
ergy loss picture is that in the medium-modified parton 
shower picture, the flow of energy from the leading par- 
ton to subleading partons is explicitly traced. Thus, such 
models make a very clear prediction that the energy lost 
from the leading shower partons leads to increased parton 
(and subsequently hadron) production at lower momenta 
in the shower. 
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7-hadron correlations [29|, [3(| HH, H2] offer a way to distin- 
guish the two scenarios. If the medium-modified parton 
shower picture is the correct way to describe the energy 
redistribution due to the medium, the additional mul- 
tiparticle production at low Pt can be experimentally 
found as a substantial enhancement of the per trigger 
yield scaled to the per trigger yield in p-p collisions (Iaa) 
in a comparatively narrow angular region (say a < 7r/5) 
back to the photon at low zt (where for each hadron 
h the relation = z^E^ links hadron energy E^ and 
photon energy £L). Here, the photon serves both as an 
indicator that a hard process has taken place and to iden- 
tify the kinematics of that process. On the other hand, if 
the energy is redistributed as a Shockwave in the medium 
[l3l [Til fla . [HI [ri| , any correlation of hadrons with the 
photon will be at thermal momenta and at large angles, 
thus it will not lead to an enhancement of Iaa when 
only an angular region of a < 7r/5 is considered to ex- 
tract the per-trigger yield of hadrons correlated with the 
photon. However, in both scenarios an Iaa below unity 
is expected at high zt since this is the kinematical re- 
gion from which energy is lost and which is also probed 
in R A a- 

In this paper, we present a systematic comparison of 
Iaa computed both in the energy loss picture and the 
medium-modified parton shower picture. For the en- 
ergy loss picture, we use the ASW quenching weights 
Note that this formalism is very well constrained 
by other data 0, [l(| EH ■ Iaa m the medium- modified 
shower picture is computed using the MC code YaJEM 
(Yet another Jet Energy-loss Model) [H, [H, [H[ using 
the ratiative energy loss (RAD) setting. Both calcula- 
tions are done using the same 3-d hydrodynamical model 
of the bulk medium evolution [33[ which is constrained 
by a large number of bulk observables. In the following, 
we consider only direct photons and assume that photons 
from fragmentation processes which would otherwise be 
a small modification to the results can be separated ex- 
perimentally, either statistically or by means of isolation 
cuts. 



II. THE MODEL 

The model has four main building blocks: 1) the primary 
hard photon and parton production, 2) the propagation 
of the away side partons through the medium and the 
parton-medium interaction 3) the evolution of the bulk 
medium and 4) the hadronization of partons emerging 
from the medium. Steps 1) and 3) are common to both 
the energy loss and the medium-modified parton shower 
picture, but 2) and 4) will be described separately for 
the energy loss picture in the ASW formalism and for 
the medium- modified parton shower picture in YaJEM. 
3) is computed in the 3-d hydrodynamical model by Bass 
and Nonaka which is described elsewhere [33|. In the 
hydrodynamical model, thermodynamical properties like 
the medium energy density e or its temperature T can be 



computed at all spacetime points. 

The framework in which all the in ingredients are merged 
is a MC code for the generation of back-to-back events 
in the medium which has been used extensively for the 
computation of hard dihadron correlations [Tol. [ill. l34l |35| 
and is described there in detail. The presentation here 
will just outline the essential physics issues. 

A. Primary parton production 

The production of a hard partons k and a photon 7 in 
leading order (LO) perturbative Quantum Choromdy- 
namics (pQCD) is described by 

dp 2 Tdyidm =2^*ifi/A(xi,Q Wj/BM )—£- 

(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k and the 
photon 7. The distribution function of a parton type i in 
A at a momentum fraction x\ and a factorization scale 
Q ~ pt is fi/A{xi,Q 2 )- The distribution functions are 
different for the free protons (36l . [37j & n d nucleons in nu- 
clei [H, [13| . The fractional momenta of the colliding par- 
tons i, j are given by xi, 2 = ^ (exp[±j/i] + exp[±y 2 ])- 
By selecting pairs k, 7 and summing over all possible 
combinations of initial partons i, j for the two contribut- 
ing channels qq — > 73 and qg — > q~f and where q stands 
for any of the quark flavours u, d, s we find the relative 
strength of the production channels as a function of px- 
In the kinematic range probed at RHIC, qg — > (77 domi- 
nates, thus quark jets are predominantly correlated with 
a photon trigger. 

For the present investigation, we require y\ = yi = 0, 
i.e. we consider only back-to-back correlations detected 
at midrapidity. In a first step, we sample Eq. |T]) summed 
over all k, 7 to generate px for the event, in the second 
step we perform a MC sampling of the decomposition of 
Eq. (JTJ according to the photon production channel to 
determine if the parton recoiling from the photon is a 
quark or a gluon. 

To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the 
nucleon (nuclear) wave function and effectively also the 
fact that hadronization is not a collinear process, we fold 
into the distribution an intrinsic transverse momentum 
kx with a Gaussian distribution, thus creating a momen- 
tum imbalance between the photon and the parton as 

PT 7 + PT k = k T 

Assuming that the photon does not interact further with 
the medium, we test at this point if the photon momen- 
tum falls into a specified trigger range. If not, a new 
event is generated, otherwise the event is processed fur- 
ther by computing the fate of the away side parton. For 
this, first the production point in the transverse plane 
must be found. 
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The probability density P(xo,yo) for finding a hard ver- 
tex at the transverse position r = {xq,uq) and impact 
parameter b is in leading order given by the product of 
the nuclear profile functions as 



P(x ,yo) = 



T A {r +h/2)T A {r -h/2) 
T AA {h) 



(2) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density p A (r,z) as T A (r) = 
J dzp A (r, z). Note that Eq. ^ may receive (presumably) 
small corrections when going beyond a leading order cal- 
culation. 

After the first step, we thus have generated a back-to- 
back event of a photon and a parton with known kine- 
matics p 7 ,Pfc and initial vertex position in the transverse 
plane ro = (xq, yo)- Partons can now be propagated from 
this position through the medium for a fixed angle with 
the reaction plane or averaged over all angles. In the 
present paper, we use the latter option. 



B. Parton-medium interaction and hadronization 
using ASW quenching weights 

For any angle <f> of the parton with the reaction plane, the 
path of the parton through the medium £(t) is specified 
by (ro, <j>) and we can compute the energy loss probability 
P(AE) path for this path. We do this by evaluating the 
line integrals 



w c (r , 



p CO poo 

= / d(<;q(0 and (ffL)(r o ,0)=/ d(q(() 
Jo Jo 

(3) 

along the path where we assume the relation 

q(C) = K ■ 2 • e 3/4 (C) (cosh p - sinh p cos /3) (4) 

between the local transport coefficient q(() (specifying 
the quenching power of the medium) , the energy density 
e and the local flow rapidity p with angle (3 between flow 
and parton trajectory. Energy density e and local flow 
rapidity p are the input from the 3-d hydrodynamical 
simulation of the medium evolution [33j . 
lu c is the characteristic gluon frequency, setting the scale 
of the energy loss probability distribution, and (qL) is a 
measure of the path-length weighted by the local quench- 
ing power. We view the parameter if as a tool to account 
for the uncertainty in the selection of a s and possible 
non-perturbative effects increasing the quenching power 
of the medium (see discussion in and adjust it such 
that pionic R AA for central Au-Au collisions is described. 
This leads to a value of K — 3.6 @. 
Using the numerical results of [Hj], we obtain 
P(AE; uj Ci R) P ath for oj c and R — 2ui^/(qL) for given par- 
ton production vertex and angle (p. We sample this dis- 
tribution to determine the actual energy loss of the away 
side parton in the event and subtract it from the away 
side parton. The away side parton k is considered com- 
pletely absorbed by the medium whenever AE > 0.9Ek- 



If the parton emerges with a finite amount of energy left, 
it is subsequently hadronized. 

Hadronization is modelled by an expansion of the frag- 
mentation function in terms of a tower of conditional 
probability densities A^(z\, . . . , z n , p) with the prob- 
ability to produce n hadrons with momentum frac- 
tions z\ , . . . z n from a parton with momentum px being 
Aj(zi, . . .Zi,pT)- The procedure is described in de- 
tail in (40| where we used PYTHIA [27| to simulate the 
shower, in the present paper we employ HERWIG [28| 
instead. We use the exact expressions for the first two 
terms and the approximations 



A 3 (zi,z 2l z 3 ,p T ) « A 2 {zi + z 2 ,z 3l p T )0(z 2 - z 3 ) (5) 



and 



A i (z 1 ,z 2l z 3l z i ,p T ) rM 2 (zi + z 2 + z 3 ,z 4 ,p T ) 
x 6(z 2 - z 3 )6(z 3 - Zi). 



(6) 



for the next two terms of the expansion. This expansion 
only becomes invalid at very low z, we will thus exclude 
this region from the discussion of results. 
Thus, for each event in the back-to-back MC code, we 
have after hadronization up to four hadrons from which 
we calculate the conditional per-trigger yield given the 
initial photon. We compute this conditional yield as a 
function of zt which is the fraction of the photon mo- 
mentum taken by each hadron. Note that due to the 
intrinsic kx smearing, zt can be above unity, unlike z. 



C. Parton-medium interaction using YaJEM 

YaJEM [H, [25| is a tool capable of computing the com- 
plete medium modified fragmentation function given a 
particular path through the medium. For the radiative 
energy loss scenario (RAD), the medium needs to be 
characterized by a parameter q which is linked to the 
thermodynamical properties of the medium by means of 
Eq. (0| , albeit with a different K — 3 which again results 
from a best fit to pionic R AA in central Au-Au collisions. 
The main physics assumption of the RAD scenario Ya- 
JEM is that the evolution equations for parton branching 
can be solved with each parton acquiring additional vir- 
tually while propagating through the medium. The gain 
of parton a in virtuality AQ^ is then given by 



AQ^= / 

Jt° 



d(q((). 



(7) 



This requires that the parton evolution which takes place 
in momentum space is linked to the medium evolution 
in position space. In YaJEM, this link is done based 
on formation time arguments. The average lifetime of a 



4 



virtual state b coming from a parent state a is computed 

as 



b = 2.4 fm 



(n) 



Ql'Ql' 



(8) 



and the actual formation time in a given event is be ob- 
tained from the probability distribution 



P{r b ) = exp 



n 
(n) 



(9) 



Based on this time information, the spatial position of 
partons in the medium is obtained by propagating them 
on an eikonal trajectory given by the shower-initiating 
parton (thus neglecting the small transverse spread of 
partons in the shower). The medium-induced virtuality 
leads to additional branching and hence an energy trans- 
fer from high pt partons into the production of multiple 
partons at lower pt- For hadronization, the complete 
parton shower is processed using the Lund model [4l| 
which is part of PYTHIA [27| . 

Due to an approximate scaling law identified in p3j , it is 
sufficient to compute the line integral 



d(q(() 



(10) 



in the medium to obtain the full medium-modified frag- 
mentation function (MMFF) Dmm(z,/i 2 ,() for a given 
eikonal path of the shower-initiating parton. Thus, for 
each event in the back-to-back MC code, YaJEM can 
provide the full MMFF, thus this is the quantity we av- 
erage after transforming it into a function of zt given the 
photon energy in the event. 



III. RESULTS 

With the two different frameworks for parton-medium 
interaction and hadronization, we are in a position to 
compute the conditional yield on the away side given 
a photon trigger both for p-p and central Au-Au colli- 
sions. To exclude the region of very low Pt, we only 
consider hadrons above 0.5 GeV. The result is shown in 
Fig.[TJ Note that there is a technical point to consider: In 
the ASW calculation, hadronization is computed inside 
the back-to-back MC simulation and thus both the sam- 
pling of geometry and energy loss and the sampling of the 
fragmentation function is done with the same number of 
events. In the YaJEM calculation, hadronization is not 
computed within the back-to-back MC code but within 
YaJEM and the full MMFF is passed into the back-to- 
back simulation. Thus, effectively the YaJEM calculation 
has a factor ^100.000 more statistics and thus vanishing 
statistical errors whereas the ASW calculation is limited 



5-7 GeV photon, ASW 
9-12 GeV photon, ASW 
5-7 GeV photon, YaJEM 
9-12 GeV photon, YaJEM 




FIG. 1: Ratio Iaa of the per trigger away side yield in cen- 
tral 200 AGeV Au-Au collisions divided by the result in p-p 
collisions given a photon trigger in the indicated range as a 
function of zt- The away side yield is extracted in a region 
in a cone with a < n/5 around the direction opposite to the 
photon. Shown are results in the energy loss picture (ASW) 
and in the medium-modified parton shower picture (YaJEM). 



in statistics at large zt — that region, however, is not so 
relevant for the present paper. 

It is immediately clear from the figure that there is a dra- 
matic difference between the expectations for Iaa for the 
two different scenarios. While the ASW scenario shows 
only a moderate rise towards low zt, the YaJEM scenario 
exhibits a dramatic increase even above unity. However, 
both scenarios agree in the high zt region as expected 
— the depletion here is in essence the same physics that 
causes single hadron Raa < 1- While there is some small 
dependence on the photon momentum, to first order the 
dramatic differences between the scenarios are indepen- 
dent of the trigger range. 

Let us first discuss the main feature, the different low 
Zt behaviour. In any medium-modified shower picture, 
there is a sum rule — if there is to be a depletion of high 
Pt hadrons in an in-medium shower as compared to a 
shower in vacuum, then energy conservation requires that 
this is compensated by an increase of low Pt hadrons. 
Thus, if the fragmentation function ratio is below unity 
in the high zt range, it must rise above unity at low zt- 
This is not so in the energy loss picture where energy 
is assumed to be distributed to extremely low momenta 
and large angles — here the ratio may stay below unity 
for any value of zt if the yield is extracted with a mo- 
mentum cut and the angular restriction on the away side 
as outlined above. 

Next, let us turn to the small rise of Iaa at low zt in the 
ASW scenario. This is in essence a kinematical effect. 
For high zt, one reqires a hadron which has approxi- 
mately (modulo the intrinsic &t) the full energy of the 
original parton. Any small amount of energy loss will 
thus immediately lead to a dramatic suppression, and 
thus one is in essence observing here hadrons from par- 



5 



tons which have not lost any energy. This is not so at 
low zj* — since the hadron energy is far from the pri- 
mary parton energy in this region, there is some room 
for energy loss, and thus one observes events in which 
there was a moderate energy loss in addition to events in 
which the parton did not lose energy, and this is reflected 
in a rise. Note that the rise is more pronounced for the 
higher trigger energy where there is even more room for 
energy loss. In any scenario in which energy loss and 
in- medium pathlength are strongly linked, as e.g. [32| . 
this kincmatical argument can be translated into a ge- 
ometrical argument leading to a transition from surface 
emission (no energy loss) to volume emission (moderate 
energy loss). However, if energy loss is a strongly fluctu- 
ating function of in-medium pathlength (as in the ASW 
scheme used here) the geometrical interpretation of the 
kinematical argument breaks down. 
Thus, there is a very clear distinction between the two 
physics pictures combined with a good insight into what 
causes the difference. 



ton. Iaa(zt) shows a clear difference between a calcula- 
tion under the assumption that energy is lost from a hard 
parton and redistributed in the medium or under the as- 
sumption that energy is transferred from hard partons 
to softer partons inside the shower during the in-medium 
shower evolution. In the case of energy loss, Iaa{zt) < 1 
is true for all zt whereas for the medium-modified parton 
shower Iaa{zt) > 1 for small zt is unavoidable. 
Preliminary data [H from both PHENIX and STAR for 
Iaa in 7-hadron correlations show no indication of a raise 
of Iaa > 1 for any zt or any trigger momentum. If these 
results are confirmed, they would constitute strong evi- 
dence that the nature of the parton-medium interaction 
is such that energy lost from an initially hard parton 
is almost completely absorbed and redistributed in the 
medium and that pictures based on an in-medium evolu- 
tion of the shower (which are in principle more complete 
than energy loss pictures) miss an essential physics pro- 
cess. 



IV. CONCLUSIONS 

We have computed Iaa the ratio of the per-trigger away 
side yield in central Au-Au collisions over p-p colli- 
sions given a photon trigger as a function zt- The 
away side yield is extracted in an angular region with 
a < 7r/5 around the direction back-to-back with the pho- 
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